Two-dimensional dissipative gap solitons 
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We introduce a model which integrates the complex Ginzburg-Landau (CGL) equation in two 
dimensions (2D) with the linear-cubic-quintic combination of loss and gain terms, self-defocusing 
nonlinearity, and a periodic potential. In this system, stable 2D dissipative gap solitons (DGSs) 
are constructed, both fundamental and vortical ones. The soliton families belong to the first finite 
bandgap of the system's linear spectrum. The solutions are obtained in a numerical form and also 
by means of an analytical approximation, which combines the variational description of the shape of 
the fundamental and vortical solitons and the balance equation for their total power. The analytical 
results agree with numerical findings. The model may be implemented as a laser medium in a 
bulk self-defocusing optical waveguide equipped with a transverse 2D grating, the predicted DGSs 
representing spatial solitons in this setting. 
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I. INTRODUCTION AND THE MODEL 

I Equations of the complex Ginzburg-Landau (CGL) type are universal asymptotic models to describe the nonlinear 
(~| ' pattern formation in dissipative media [Ij]. They also find direct (rather than asymptotically derived) realizations 
' ^ in nonlinear optics as models of laser cavities [Ij. Objects of fundamental interest predicted by the CGL equations 
^ ,. are solitary pulses (SPs), alias dissipative solitons. They represent, in particular, temporal pulses generated by fiber 
' lasers 0,0,131 and, in an altogether different physical context, patches of traveling- wave thermal convection in narrow 

'. channels @. 

^ ' The simplest CGL equation is based on the cubic nonlinearity. Exact one-dimensional (ID) SP solutions are known 
\ in that case [6|, but they are unstable, as the respective equation includes the linear gain, compensating the cubic 
C i loss, which destabilizes the zero background around the pulses and leads to the formation of a chaotic "gas" of SPs 
' 0- ^ well-known modified equation that can support stable SPs includes a combination of linear and quintic loss 
, terms and cubic gain, therefore it is called the cubic-quintic (CQ) CGL equation. This equation was introduced (in 
OO the 2D form) in Ref. Q. Stable SP solutions to the ID version of the CQ-CGL equation were first predicted by 
means of an analytical approximation based on the power-balance analysis for solitons of the nonlinear Schrodinger 
' (NLS) equation Then, these solutions were explored in detail by means of numerical methods [l^]- More general 
models, such as linearly coupled systems of CQ-CGL equations [11] and the complex Swift-Hohenberg equation with 
the CQ nonlinearity |12l |. were introduced too. 

In addition to the solitons of the NLS type, a generic species of solitary waves in conservative media is represented 
by gap solitons. They are well known in nonlinear optics, as temporal solitons in fiber Bragg gratings, and in 
Bose-Einstein condensates (BECs) with repulsion between atoms, which are loaded into in optical-lattice trapping 
potentials. A fundamental feature, from which the name of the gap soliton derives, is that its wavenumber (in terms 
of optical models) must belong to a finite bandgap induced by the effective periodic potential. In the context of both 
the fiber gratings [l3l[T^ . [T5l and matter waves (BEG) [H, [l3|) gap solitons were predicted theoretically, including 2D 
gap solitons in BEG [l^ [l^ , and 2D solitons with embedded vorticity |2(|. The creation of optical and matter- wave 
gap solitons was reported in the ID geometry - in short fiber gratings [22l |. and in the condensate loaded into an 
optical lattice combined with a strong transverse trap [l^ , respectively. 

A generalization of the above concepts, aiming to predict dissipative gap solitons (DGSs), was proposed recently 
[23 | . The respective version of the CGL equation combines a periodic potential and the set of the CQ loss and gain 
terms: 



dt " ^29^ 

*e(-7i +721^1' "731^-1")^. (1) 
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Here, ip may be considered as the local amplitude of the electromagnetic wave (in the case of an optical model), A 
and n/qo are the strength and period of the potential {x = L/2 is the midpoint of the system), and coefficients €71, 
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£73 cind £72 represent the linear and quintic dissipation and cubic gain, respectively. The conservative cubic term in 
Eq. ^ is defined with the self-defocusing sign, which is relevant to gap solitons. 

A physical realization of Eq. ([1]) pertains to a planar self-defocusing optical waveguide, in which the periodic 
potential may be induced by a transverse grating (periodic modulation of the refractive index). In particular, a 
self-defocusing nonlinearity, in the combination with photoinduccd transverse lattices, can be implemented in pho- 
torefractive crystals [1^. As concerns the CQ loss/gain terms, they may actually represent a combination of the 
linear amplification and saturable absorption, which is a common setting in laser cavities 0, l26| . In terms of this 
interpretation, variable t in Eq. ([1]) designates the propagation distance, while x is the transverse coordinate [hence 
the SP solutions to be produced by Eq. ([!]) will be spatial solitons [2J|]. It is also relevant to mention the model 
combining the standard fiber-Bragg-grading part and the CQ combination of dissipative terms, in which temporal 
DGSs were investigated [131, including interactions between them (28j . 

The analysis reported in Ref. j24| revealed, by means of approximate analytical and direct numerical methods, 
the existence of three families of stable DGSs in the first finite bandgap of the respective linear spectrum: loosely 
and tightly bound static solitons, and a family of breathers between them. All the families were found close to the 
border between the finite bandgap and Bloch band separating it from the semi-infinite gap, the tightly and loosely 
bound DGSs being located in spectral regions where the Bloch band is, respectively, very narrow or relatively wide. 
Stable dark solitons were also found in the model, and the mobility of dark and loosely-bound bright solitons was 
demonstrated in it (collisions between mobile solitons are quasi-elastic). 

The objective of the present work is to introduce fundamental and vortical DGSs in two dimensions. To this end, 
we consider the 2D extension of Eq. ([1]), 
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This equation may be interpreted as governing the propagation of electromagnetic waves in a bulk medium with the 
self-defocusing nonlinearity and other ingredients included in Eq. ([T]), assuming that the potential periodic in x and 
y is induced by the 2D transverse grating, with t again having the meaning of the propagation distance. In Section 
II, solutions for fundamental gap solitons are found, in parallel, in a numerical form [as attractors of Eq. Q], and 
by means of an analytical approach, which combines a variational approximation (VA) for the shape of the solitons 
and the balance equation for their total power. In terms of the corresponding linearized equation, the soliton family 
belongs to the first finite bandgap. Further, in Section III we report the existence of stable vortex solitons (with 
topological charge 1), built as complexes of four peaks, with the phase shift of tt/2 between adjacent ones, and an 
empty site in the center {rhombus- shaped vortices [29| . alias on-site ones). The vortex solitons are constructed as 
stable numerical solutions, and are also obtained by means of the analytical approximation. Thus, the results reported 
in this work predict the existence of stable spatial solitons, both fundamental and vortical ones, in laser cavities based 
on self-defocusing bulk media with transverse gratings. 



II. FUNDAMENTAL SOLITONS 



Stable solutions to Eq. ^ in the form of 2D fundamental (zero-vorticity) solitons can be readily found, in the first 
finite bandgap, as attractors of the CGL equation, by dint of direct simulations starting with an appropriate initial 
configuration. Typical examples of such tightly and loosely localized gap solitons are displayed in Fig. [TJ The 2D 
solitons become looser with the decrease of the lattice period, tt/qo, as the soliton spreads over a larger number of 
lattice cells. 

The shape of the fundamental solitons is determined by the equation for complex function (j){x,y), obtained by 
the substitution of 'ip{x,y,t) ~ e~'^*(/) (a:, y) in Eq. ([2]), where — ^ is the soliton's propagation constant, in terms of 
the optical model. To approximate the soliton's shape in an analytical form, we start with the equation for (p {x, y) 
without the dissipative terms {e = 0). The latter equation can be derived from Lagrangian A = J J Cdxdy, with 
density 

} i^r- (3) 

In this approximation, (p(x, y) may be assumed real, therefore, following the lines of Refs. and [3l|, we adopt the 
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FIG. 1: Numerically generated profiles of stable fundamental tightly and loosely bound 2D gap solitons, \i};{x,y)\, found, 
respectively, at go = 1 (a), and at qo — 1.75 (b). Other coefficients are 71 — 0.5,72 = 2,73 = 1, e = 0.05, and A — 2. 



variational ansatz as 
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where aniplitude B and width a ^ are free parameters, the norm of ansatz (jj]) being 
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(in terms of the laser-cavity model, N is proportional to the total power of the light beam). 

Ansatz implies that constant q may be different from go in Eq- ©■ In fact, for the description of tightly-bound 
(strongly localized) solitons it will be sufficient to fix g = go, which is not surprising, as the mismatch between the 
periodic functions in the ansatz and the periodicity of the potential, accounted for by g 7^ go, is not crucially important 
for the well-localized waveforms. On the other hand, applying the VA to loosely bound (weakly localized) solitons, 
we will treat g not as a variational parameter, but rather as a "phenomenological" fitting constant. We also tried 
an extended version of the VA which subjected g to variation too, but it yielded essentially less accurate results. 
This conclusion may be explained by the difficulty in approximating the complex shape of loosely-bound solitons 
by a simple ansatz. Similar situation are known in other applications of VA to the description of complicated wave 
patterns, where some parameters should be still treated as variational ones, while others are reserved for direct fitting, 



cf., e.g., Ref. [32 1 



The substitution of ansatz ^ into Lagrangian density ([3]) and the integration yield the effective Lagrangian, 
expressed in terms of A^, a and g: 



= -9 /a 
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Variational equation dA/da = for fixed g determines a as a function of N (as said above, g is not a variational 
parameter, hence we do not add equation dA/dq — 0). 

At the next stage of the analysis, we restore the dissipative terms in Eq. ^ and, treating them as small perturba- 
tions, derive a straightforward evolution equation for the total norm. 



~dt 



(7) 



4 



(a) 



(b) 



(c) 



(d) 




1.5- 




2- 


■3- 1- 




1.5- 
" 1- 


0.5- 




0.5- 


0- 




0- 



-5 5 

X-U2 




FIG. 2: (a) Numerically found amplitude B of tightly bound dissipative gap solitons (diamonds) and the respective analytical 
prediction (lines), versus cubic gain 72. Other coefficients are A = 2,qo = l,e = 0.05, 71 = 0.5,73 = 1, and q — qo = 1- 
(b) The comparison of the numerically found and analytically predicted cross sections of the soliton profiles, \ip{x,y — L/2)\, 
for q — qo — 1, e = 0.05, 71 — 0.5,73 — 1, and 72 = 2 (the case of relatively tightly-bound profiles), (c) Numerically found 
amplitude B of loosely bound dissipative gap solitons and the respective analytical prediction, obtained with fitting parameter 
q = qo = 1.75, or q — 0.6go = 1.05, versus cubic gain 72, for qo = 1.75. (d) The comparison of the numerically found and 
analytically predicted cross sections of loosely bound soliton profiles for go = 1.75, with the choice of q — 1.05. 



where the following coefficients were calculated as per ansatz @: 
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According to Eq. the balance condition for the norm in the stationary state, dN/dt — 0, yields relation 

727V2 =7iA^ + 73A^3- (10) 

Figure l^a) displays the amphtude of the DGS, B for — 1, versus cubic gain 72, as predicted by the analytical 
approximation with q = qo, i.e., obtained from a numerical solution of equations ([5]) and (jlO[) . along with the same 
dependence generated by direct simulations. It is seen that the analytical approximation provides good overall 
accuracy for the solitons which feature a (relatively) tightly bound shape (therefore, it is sufficient to set q = qo, i.e., 
one does not need an extra fitting parameter in this case), although the approximation predicts the existence of the 
solitons at 72 < 2.15, while the direct simulations yield stable DGSs at 72 < 2.7. For 72 > 2.7, the simulations produce 
delocalized solutions, as the loss terms cannot compensate the gain in that region. The trend to underestimating the 
existence area for the DGS by the analytical approximation is generic, persisting throughout the parameter space. 

For the same case of g = go = 1, Fig. Wi^) displays comparison of the soliton's cross-section profile, \ip{x, y = L/2)\, 
and the respective analytical approximation. For given values of the parameters, the latter one is \il){x,y = L/2)\ = 
1.53e-°-4"a; cos(a;-L/2), seeEq. Q. It is seen that the analytically predicted and numerically found profiles overlap 
well near the soliton's center, but points at which (t>{x) crosses zero are shifted. 

Figure [^Jc) displays the amplitude of loosely bound solitons versus 72, for qo = 1.75. The plot shows both the 
straightforward prediction of the analytical approximation, obtained with q = qo — 1.75, and its "phcnomenologically" 
adjusted modification, provided by choosing q — O-Gqo = 1.05. In the direct numerical simulations, the respective DGSs 
exist for 1.95 < 72 < 2.95, while the analytical solution with q = qo = 1-75 exists in the region of 1.64 < 72 < 1.75, 
which does not overlap at all with its numerical counterpart. However, choosing the fitting parameter to be g = 1.05 
allows us to fit the region of the existence of the analytical solutions to the numerical one fairly well. Accordingly, 
the comparison of the numerically found soliton's cross-section profile, \4>{x,y = L/2)\, and the respective analytical 
approximation with q ~ 1.05, which is \tp{x, y = L/2)\ = 1.46e~'^'^^^^ cos 1.05 (a; — i/2) in the present case, is displayed 
in Fig. Hd). 

Figure [3] shows another set of global characteristics of the DGS family, viz., the amplitude and propagation constant 
(B and /i) versus the grating's wavenumber, as found in the numerical and analytical forms. In particular, panel[3fb) 
clearly shows that the entire DGS family falls into the first finite bandgap of the underlying linear spectrum. The figure 
also demonstrates the same trend as mentioned above, namely, that the analytical approximation, while predicting 
generally correct characteristics of the soliton family, underestimates its existence range: it does not yield solutions for 
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FIG. 3: Amplitude B (a) and propagation constant /i (b) of fundamental solitons versus the grating's wavenumber, 2go, for 
A — 2 and e = 0.05, 71 = 0.5,72 = 2, 73 = 1. Diamonds and bold dashed curves show, respectively, direct numerical results 
and predictions of the straightforward analytical approximation, obtained with q = qo (no fitting). The thin dashed curve in 
(b) is the lower border of the first finite bandgap [as found with e = 0, i.e., in the conservative version of Eq. ([2j], to which the 
soliton family belongs. 
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FIG. 4: Stable fundamental gap solitons exist above the continuous curve in the (71, A) plane. Other parameters are 72 = 
2, 73 = 1, e = 0.05, and go ~ 1. 



go > 1-2, while the direct simulations converge to DGSs up to g = 1.85 (numerically found solutions get delocalized 
at go > 1-85). 

The existence range for the stable DGSs in the plane of parameters which control the linear part of the model, viz., 
loss coefficient 71 and grating strength A, as found from the numerical results, is displayed in Fig. [H To the right of 
the border shown in Fig. |4](i.e., if the linear attenuation is too strong), numerical solutions decay to zero, while to 
the left of the border (if the attenuation is too weak) the solution undergoes a delocalization transition. It is relevant 
to mention that, unlike DGSs in the ID counterpart of the present model [24*1, no 2D soliton was found to be mobile 
[the application of the kick, i.e., multiplication by exp (iKx), fails to set any 2D soliton in persistent motion]. 



III. SOLITARY VORTICES 

We have also found dissipative gap vortex solitons with topological charge 1. To generate the vortex soliton, we 
have performed numerical simulation of Eq. (2) starting with the initial configuration 



■)pix,y) = B 



X exp ■ 



^"2 



g X 
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L 



q y 



(11) 



where the pre-exponential multiplier accounts for topological charge 1 of the vortex. Typical examples of stable 
vortices are displayed in Fig. [51 at the same values of parameters for which stable fundamental solitons were shown 
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FIG. 5; Profiles of vortex gap solitons, \tp{x,y)\, produced by tfie direct simulations at go = 1 (a) and 1.6 (b) for 71 — 0.5,72 = 
2,73 = 1, e = 0.05, and A = 2. 
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FIG. 6: (a) Tiie same as in Fig. [2ja) (at the same values of parameters), but for the amplitude of the rhombus-shaped vortex, 
(b) The same as in Fig. EJb), but for the cross-section profiles of the vortex, (c) The numerically found vortex' amplitude for 
qo — 1.6, and the respective analytical approximation (the dashed line) adjusted by choosing q = 0.7qo = 1.12. (d) The same 
as in (b), but for qo = 1.6 and q = 1.12. 

in Fig. [1] Figure ini^a) displays a rhombus- type vortex soliton for qq = 1, which is composed of four peaks with the 
phase shift of tt/2 between adjacent ones. Figure [S^b) displays a more loosely bound vortex soliton, for qq = 1.6. 
The solitary vortices feature a transition from the tightly bound shape to the loose one with the decrease of the 
period of the underlying grating, tt/qo- In any case, main peaks which form vortices are located, approximately, at 
{x,y) = {L/2 + Tr/qo,L/2),{L/2- n/qo,L/2),{L/2,L/2- n/qo) and {1/2,1/2 + n/qo). The vortex slowly becomes 
delocalized at still larger values of go, for instance at go = 1.75, if other parameters are fixed as in Fig. [5] 

Getting back to the example of a tightly bound vortex obtained at go = 1 , which is built of the four peaks which are 
well localized around four potential minima, (x, y) = (L/2-f tt, £/2), (L/2 — tt, i/2), (L/2, L/2 + -k) and {L/2, L/2 — tt), 
one may approximate each peak by ansatz ^ with q = qo. Then, the entire vortex may be approximated as 

{x, y) — B cos{x — L/2) cos(j/ — L/2) x 
{exp[-a/2{(a; - L/2 - n)^ + {y - L/2)^}] + i exp[-a/2{(a; - L/2)^ + {y-L/2- Tr)^}] 
- exp[-a/2{(a; - L/2 + Tr)^ + (y - L/2)^}] - i exp[-a/2{(a; - L/2)'^ + {y-L/2 + n)^]} (12) 

where B and a are taken as predicted by the analytical approximation for the fundamental soliton. 

Figure [6lJ a) shows the dependence of the vortex' amplitude on the cubic-gain coefficient, 72, as found from the 
direct simulations for go = 1, i.e., for the family of relatively tightly-bound vortices. The dashed curve in Fig. [5]^ a) 
represents the corresponding analytical approximation for B in expression (|12p. which is found to agree with the 
numerical results. Further, Fig. [5] (b) displays a typical example of the numerically and analytically found cross- 
section profiles of the tightly-bound vortex, which also demonstrates good agreement between both [the dashed curve 
is li^ix)] = 1.532| cos(a; - L/2)[exp{-0.413(a; - L/2 - n)^} - exp{-0.413(a: - L/2 + 7r)2}]|, as per Eq. P^]. 

Proceeding to loosely bound vortices. Fig. [6Kc) shows the dependence of the vortex' amplitude 
on 72 for go = 1.6. The dashed curve in the same panel represents the respective analytical ap- 
proximation for B, adjusted by choosing g = 0.7go = 1.12 in ansatz Further, Fig. Eld) dis- 
plays a typical example of the numerically and analytically found cross-section profiles of the vortex. 



7 



also for qo = 1.6 and q = O.lqo = 1.12. The dashed (analytical) profile corresponds to = 
B |cos [q{x - L/2 - n/qo)] exp [-(a/2) (a; - L/2 - n/qo)^] - cos [q{x - L/2 + n/qo)] exp [-(a/2)(x - L/2 + 7r/go)^] |, 
where B — 1.461 and a = 0.339 are values predicted by the analytical approximation for the fundamental DGS. 
Good agreement between both profiles is observed. 



IV. CONCLUSION 

In this work, we have introduced a model combining the CGL equation in two dimensions with the CQ (cubic- 
quintic) combination of gain and loss terms, self-defocusing nonlinearity, and the 2D periodic potential. The model 
gives rise to stable DGSs (dissipative gap solitons), both fundamental and vortical ones. The DGS families were 
found in numerical and approximate analytical forms in the first finite bandgap of the model's linear spectrum. The 
analytical approximation, based on the VA (variational approximation) for the soliton's shape and balance equation 
for the total power, yields results which turn out to be quite accurate in comparison with the numerical findings (for 
the fundamental and vortical solitons alike) , except for the fact that the analytical approximation underestimates the 
region of the existence of stable DGSs. The model considered in this work can be realized as a laser cavity in a bulk 
self-defocusing optical waveguide equipped with a transverse grating. 

The analysis of the model can be developed in other directions. In particular, it may be interesting to construct 
solitary vortices with multiple values of the topological charge (stable 2D gap solitons with embedded vorticity 2 were 
found in the conservative model 21]). A challenging problem is to find DGSs and vortices in higher finite bandgaps 
of the underlying linear spectrum. 
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